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Abstract 

We compute the Yang-Mills vacuum wave functional in three dimensions at weak coupling with 
0(e 2 ) precision. We use two different methods to solve the Schroedinger functional equation. One 
of them generalizes to 0(e 2 ) the method followed by Hatfield at 0(e) The other uses the weak 
coupling version of the gauge invariant formulation of the Schroedinger equation and the ground 
state wave functional followed by Karabali, Nair, and Yelnikov [2J. We compare both results and 
discuss the differences between them. 

PACS numbers: ll.10.Ef, ll.10.Kk, 12.38.-t 



I. INTRODUCTION 



The determination of the ground-state (or vacuum) wave functional of QCD, *&[A], is 
tantamount to solving QCD, as any observable (for instance the static potential or the 
spectrum of the theory) can then be obtained by the computation of the expectation value 
of the appropriate operator. Even if the exact solution is not known, properly chosen trial 
functions may give valuable information of the vacuum using variational methods (see for 
instance [3]). 

We are still far from obtaining the exact ground-state wave functional of QCD. Even 
obtaining approximated expressions is very complicated. This is also true in the weak cou- 
pling limit. One reason follows from demanding the wave functional, besides of satisfying 
the Schroedinger equation, to be gauge invariant. This constraint is imposed by the Gauss 
law, and makes that standard quantum mechanics perturbation theory can not be applied. 
A procedure to overcome this problem was devised in the case of SU(2), and to 0(e) in the 
weak coupling expansion, in Ref. pQ. An alternative procedure has also been considered in 
Ref. [3] and worked out to the same order in e. 

In this paper we are interested in the three dimensional version of QCD without light 
fermions (i.e. Yang-Mills theory or gluodynamics) . The method outlined by Hatfield pQ 
can also be applied to the three dimensional case and a general group SU(N) without major 
modifications. We do so in sec. [n]and obtain the 0(e) expression for a general group SU(N) 
in three dimensions. The result is in perfect agreement with the expression obtained by 
the expected replacements of the four dimensional result obtained in Ref. pQ. The results 
obtained with this method satisfy the Schroedinger equation by construction but not nec- 
essarily the Gauss law, though it can be explicitly shown that it does at 0(e). We then 
compute the 0(e 2 ) wave functional in what is a completely new result. Again, this result 
satisfies the Schroedinger equation by construction but at this order it's not possible to 
explicitly check the Gauss law, due to the complexity of the resulting expressions. The re- 
sulting wave functional is explicitly real (as expected for the ground-state functional) and 
we name it ^g^LA]. 

The fact that gauge invariance can not be guaranteed in general is one important drawback 
of the previous method. The reason is that the Gauss law is only implemented partially for 
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some terms in some intermediate expressions. Moreover, even this partial implementation 
of the Gauss law is difficult to automatize, as at each order it has to be tailored somewhat. 

One solution to the previous problem would be to reformulate the Schroedinger equation 
in terms of gauge invariant variables. One such formulation was originally worked out in Refs. 
[HHH] (for some introductory notes see [10J ) and, more recently, in Ref. [2J, where a modified 
approximation scheme was devised. The authors use a change of field variables, which become 
complex, to simplify the problem. Even though the original motivation of those works was 
to understand the strong coupling limit (the opposite limit we are considering in this paper) 
and confinement in three dimensions, it is not difficult to see that the approximation scheme 
worked out in Ref. [2J could be easily reformulated to provide with a systematic expansion 



of the weak coupling limit. We do so in sec. Ill of this paper and compute the ground-state 
wave functional to 0(e 2 ). The vacuum wave functional is a function of the gauge invariant 
variables J a , which we then transform to the original gauge variables A a . The resulting 
expression is gauge invariant by construction and also satisfies the Schroedinger equation by 
construction. We name it \l/ar[A] = ^gi [J {A)]- Nevertheless, the explicit expression has the 
very unpleasant feature of having a non-trivial imaginary term. 

We have then obtained two different expressions for the wave functionals: \&gl[A] and 



\&(3jL4], which actually look completely different. We compare them in sec. IV. At O(e) 
it is possible to show after several manipulations and using the symmetries of the integrals 
that both of them are equal (so at this order they are real and gauge invariant). Such brute 
force approach happens to be unfeasible at 0(e 2 ) due to the complexity of the expressions. 
We need an organizing principle for the comparison. The approach we follow is to rewrite 
\&gi<[-A] in terms of the gauge invariant variable J a and a gauge dependent field 9 a . All 6 a 
dependent terms should vanish if ^^[A] is going to satisfy the Gauss law and we explicitly 
show that this happens. This means that both *&g?l[A] an d ^g/[AI are gauge invariant. We 
would then say that both should be equal, since both satisfy the Schroedinger equation. We 
actually find (after a rather lengthy computation) that both are almost but not completely 
equal. The difference is proportional to a bilinear real term. This is puzzling but there is 
a reason behind it: \&gl[A] and ^/[A] satisfy "different" Schroedinger equations. ^g^LA] 
was obtained using the unregulated Schroedinger equation, whereas ^/[A] was obtained 
after the Schroedinger equation in terms of J a variables was regularized. In this last case, 



regularization produces an extra term in the Schroedinger equation, producing in turn an 
extra term in the wave functional. Irrespectively of the above this comparison allows to 
rewrite ^^[A] in an explicitly real form. This is by far non-trivial, as the initial ^^[J] 
was explicitly complex and dependent on complex variables. In particular there is a delicate 
cancellation between terms such that, after transforming this expression back to real vari- 
ables, the wave function becomes real (actually in our comparison we work the other way 
around and transform ^^[A], which is real, in terms of the complex variables). This is an 
important test of several parts of the computation done in Ref. [2]. 

We believe that the weak coupling reformulation of the approach followed in Ref. [2] 
can be helpful to understand the meaning of the partial resummations performed in the 
approximation scheme used in this reference, though we do not explore this issue in this 
paper. Our 0(e) or 0(e 2 ) wave functional can also be used to test different trial functionals 
that one can find in the literature claiming to have the proper weak and strong coupling 
limit. Typically they reproduce the leading order weak coupling expansion but not the 0(e) 
corrections. This is certainly the case with covariantization approaches where the exponent 
of the wave functional is approximated by a bilinear term in the B fields (see for instance 
[TT| fT2]). Therefore, our results can hint to how those trial functions could be improved to 
correctly incorporate corrections in the weak coupling limit. 



II. DETERMINATION OF ^ GL [A] 
The Yang-Mills Lagrangian reads 



where 



C = - 1 -G^ a G% 



G% = d,A a v - d u Al + et hc A\Al , (2) 



eGfj, v = [D^,D U ], D„ = <9 M + eA^ A„ = —iT a A a ^ G^ = -iT a G% u , T a = \ a /2, \ a are the 
Gell-Mann matrices, and [T a , T b ] = if abc T c . 

We will work in the Hamiltonian formalism and partially fix the gauge to A = 0. There- 



fore we work with the A = (Ax, A 2 ) components only and 



D = V + eA, (3) 

B a = l -e ]k (d 3 A k - d k A 3 + e[Aj, A k ]) a = VxI«+ 6 -f ahc A b x A c , (4) 

where A x B = tijAiBj, V« = di = d/dx % (for simplicity, we use the metric j]^ = 
diag(— 1, +1, +1), so there is no sign difference between upper and lower spatial indices), 
and B = -%T a B a . 

In Ref. pQ the wave functional was computed to 0(e) at weak coupling. It is possible to 
generalize the method used in this reference. We do so here and compute the ground state 
wave function to 0(e 2 ). The ground state wave functional has to satisfy the Schroedinger 



equatior^ 1 



\ I \-JT7^ • JTT^ + Ba ^ Ba W ) ^l[A] = E* GL [A] , (5) 



and the Gauss law constraint 



( D ■ LT*ol[A] { V • + ef abc A b ■ -* ) *« L [A] = . (G) 

SA n 5A C 



Because we are talking of the ground state we expect the wave function to be real and to 
have zero nodes. Therefore it can be written as the exponential of a functional F[A], which 
does not diverge for finite A: 



m GL [A] = e - FcL[A] = e -F { GM-eF£i[A]-e*Fgl[A]+0(e*) _ ^ 



This particularly simplifies the Gauss law, which then can be written as 

^ • TT + eflbcAb ' TT ) F ol[A] = . 
oA a 6A C J 



1 In the following we use the notation (d = 2): J x = J d d x, J^ = J j^Jjd, ${k) = (2n) d 8^ d \k) , and so on. 

2 Note that the ground state energy can be normalized to zero by moving it to the left-hand side of the 
equation and absorbing it in the B 2 term as a counterterm. 
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A. Order e' 



FqI can be obtained in several ways. It is equivalent to solving the Schroedinger equation 
of the free theory with the free Gauss law, in other words, N 2 — 1 replicas of QED without 
light fermions. In order to solve these equations it is convenient to rewrite them in momentum 
space using: 

A(x) = I A(k)e^, = [ -j^re**. (9) 

We then have for the free-field Schroedinger equation 

\ [ \-7jL\ ■ VY^IS + (* x x Ia (-V)) = EV>*®[A] , (10) 

2 Jf \ 5A a (k) 5A a (-k) J 

giving the following equation for the leading order term of the wave functional exponent 

/ itFtk ■ vfr^ = A* x Aa ^ k x Aa (~^ ■ ( u ) 

J$ 5A a (k) 5A a (-k) Ji 
The free-field Gauss law reads: 

k ■ g L L = 0. 12 
SA a (k) 



Eq. (J 1 1 [) suggests FqI to be quadratic in A: 



A«(k)A"(-k) 9ij (k) 



(13) 



The tensor structure of gij{k) can be fixed by the free-field Gauss law, Eq. (12), which 
demands that only the transverse component should show up. Therefore 



9ij(k) = g(k)Vij(k) , 



(14) 



where Vy = 5{j — kikk/k 2 is the projector to the transverse component. We can now solve 



Eq. (11) and determine g(k). As the equation is quadratic there are two solutions, of which 



we take the one that leads to a normalizable wave functional, which is 

F [ gI[A] = \ f^r(kx A a (k))(k x A\-k)) , (15) 

where Ek = \k\. A detailed explanation of this derivation can be found in sec. 10.2 (see also 
sec. 11.2) of [13]. One can see how the implementation of the Gauss law is not trivial even 
in the free-field case. 



B. Order e 



At 0{e) the Schroedinger equation splits into two equations (organized by powers of A): 



5F§l[A] 6Fg>[A] _i 
1 5A a (-k) 5A a (k) 2 



f abC [ $ fe^l & x ^ a (£i))(^(£ 2 ) x A c (h)) , (16) 



S 2 F^ L [A] 



If 5A a (-k)5A a (k) 
and the Gauss law constraint reads^] 



5A a (k) 



k ■ 



-if abc 
-if abc 



ft Pi • 

5A a (p 2 ) 



P/l,P/2 
1 

w 



P2 + k) 



(p x A^i-k-p)) (p x A c {p)\ . 



Using Eq.(15) the left-hand-side of Eq. (16) can be rewritten as follows: 

8Fg:l[A]^ 

6Fgl[A] 



^j(px A a {p)) [px 
\P\ V 5A-(p) 



1 yiA'M 5Fg[[a] 

\P\Y [ A{P) 5AHp 



p. A a ( 



V ■ 



5A a (p) 



(17) 



(19) 



where the second term of the right-hand-side is known because of the Gauss law. 

We are now in the position to obtain F^: We profit from the fact that the kernel can be 



Note that in d = 2: A c (-k - p) ■ (px (px A b (p) 



px A b {p)\ {px A c {-k-p)). 



Other useful relations are (k ■ A)(k x B) — (k x A)(k ■ B) = k 2 (A x B) and eyey = S^Sji ~ SuSjk- 
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taken to be completely symmetricn under the interchange of any two fields Ai a . x ., A ia . x .. 

Th e Kfore , the L ty of \ M . can be related L the densit ;;; ^ 

More specifically, if for a functional F A ai (ki), . . . , A an (k n ) of n fields we have 



\p\ A a (p) 



SF[A) \ 
5A a (p) I 



D 



A ai (kx), . . .,A an (k n ) 



(20) 



then 



F[A] 



1 



-D 



+ • • • + \k r , 
With this we finally obtain 



A a ^h),...,A a -(k n ) 



(21) 



Fgi[A} = ir be l ^(££ 

JlyflMM \ i=1 

1 



-(h x A»(A; 1 ))(i b (A; 2 ) x A c (k 3 )) 



(E?N)NN 



2(£ l 3 N) 

(*i • I a (£i))(£ 3 x j4*(£ 2 ))(£ 3 x A c (k 3 )) 



(22) 



which is the three dimensional version of Hatfield's result (except for a different sign con- 
vention for e). 



C. Order e 2 

At 0(e 2 ) the Schroedinger equation leads to the following equality 
1 r /A2p( 2 ) afW ip(0) AP (2) -, \ 

1 / I GL GL GL o or GL GL , 1 rabcrade/ 7b v ac\/ Td ^ Te\ \ _ n / 0Q \ 

2Al v (^~(^(^~ 2 (^(^ + 4 / ; (^x^)(^x^)J-0- (23) 

At this order .F^ can have contributions from four, two and zero fields (there are no contri- 
butions with three or one field): Fq\ = FqI Fq£\ There is no need to compute 
FqJj, as it just changes the normalization of the state, which we do not fix, or alternatively 



can be absorbed in a redefinition of the ground-state energy. Then, Eq. (23) can be split 



4 Any term antisymmetric in any of the two indices will vanish when multiplied by the gauge fields. This 
means that the kernel is not completely determined, as such terms can always be added. 
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into the term with two and four fields: 




(0) X 17(2,4) 



6F$> SF 



SF§1 SF§1 

(8 At) (8 At) ' "{8 At) (5 A 



GL 



£r + -f abc f ade (A b x A c )(A d x A e 



(24) 



and 




5 2 F 



(8Aty (8At) (SAf) 



0. 



(25) 



FqI and Fq[ have already been determined (see Eqs. (13) and (22)) and can be inserted 



into Eqs. (24) and (25), but we still have to implement the Gauss law, which at this order 
reads 



fp(2,4) , r P (l) 



SA a (k) 



8At(p 2 



P2 + k) 



(26) 



5A a (k) 



0. 



(27) 



(2 4) (2 2) 

One first solves Eq. (24) and determines Fq£ . Afterwards Fq£ is fixed by Eq. (25). In 



order to obtain Fq£* the procedure is similar to the one used for FqI- The dependence on 



(i) 



Fq^ is encoded in the 2nd term of Eq. (24), which we rewrite in the following way 



8F%£\A] 



— (pxA a (p)) px 
IP] \ 8A a 




5A a {p) 



p-A a 




(2f 



Once again the second term on the right-hand-side is given by the Gauss law, which allows 



us to isolate F. 



(2,4) 



GL 



As above we use the fact that the kernel can be taken to be com- 



pletely symmetric under the interchange of fields Ai i(HjX ., which lets us relate the density of 

L\p\ (Hp) 



,(2,4) 
GL 



SA a (p) 



with the density of Fq^ [A] and we finally obtain 



7 (2,4) 
GL 



SF 



(i) 



GL 



2 J***** EJdfel + 1*1) \ 5A ?® 



[fci,fc 2 



SF, 



(i) 



GL 



SAt(-p) 



[91,92] 



cbib 2 c 



Hqi + Q2-P) 

# + Wi\)Wl\ 



Xb-2 



SF W -. - 
SA c {p) 



9 



8 JthMofi Ei(N + Iftl) 



k E; < +" ^1 x ^a 2 (^)) (ft) X ^ (g 2 )) , (29) 



which explicitly reads 



tti(2,4) r-abc rede 
GL = 



f abC f Cde [ $ + $ ) , r , , r / 

J %,Hi,q/L,qfl \ i J \Kl\ + |K2| 4 



1 



2(|fci| + N + |fci + fc 2 |)(|gi| + |g 2 | + |gi + ? 2 |) 



C2| + |?l| + 1 92 

^ d (gi) x A e (q 2 ) 



-^Aix + Aial^fci) x^ 2 ) 



+ |fcl ,j", fc2| (fci + k 2 ) x l a (£i)(£ 2 • ^(jfe 2 )) + ^±M-^ (jb! ■ l a (^))(fc 2 x A 6 (/c 2 )) 



1*2 

x x (>(£ 2 ) • A c (g 2 )) 



Ml** 



1 



2^ 2 • l e (g 2 ) - ^ig 2 • A e (g 2 ) 



l9ip 2 



(fci • A a (h))(k 2 x ^(A; 2 ))(gi x 



1 



+^(k 1 -A a (k 1 ))(k 1 + k 2 )xA b (k 2 ) 
\ki\ Iml 

+T*^itfi x A d {&)){& + 92) • A\q 2 ) 
2 (gi + gg) • qi 

2ki x fc 2 



(91 + 92) x A d {q x ){q 2 - A e {q 2 )) 



\ki\\k 2 \\qi + q 2 \\q 2 \ 



(h ■ A^h^h + k 2 ) x tffyiqi x ^ d (gi))(g 2 • ^ e (g 2 )) 
(fci • l a (£i))(£ 2 x ^(fc 2 ))(£ x ^ d (gi))(g 2 x l e (g 2 )) 



(h x A a (A; 1 ))(A; 1 + k 2 ) x A fe (A; 2 )(g 2 x A d (gi))(g 2 x A e (q 2 )) 



9i + 92II92I 

; (fc 2 x l a (^))(fc 2 x #(£ 2 ))(<f 2 x A d (g 1 ))(92 x l e (g 2 )) 



I&2II92I 



1 



+ ^ (X"^) x #(£ 2 )) x X c (g 2 ) 



- — ^— — - — (fci " ^ a (*i))|§(*i + fe) x 1^ 2 ) (l d (gl) x A\q 2 )) 
-(gl x A d (ft) (>(£ 2 ) x l e (g 2 )) 



1 



(fci + fc 2 ) x A b (/c 2 )(gi + q 2 ) x A'Vi)^ • A e (q 2 )) 



I91 + 92II92I 



(q 2 x i*(Al 2 ))(gi • A d {q x )){q 2 x X e (g 2 )) 



I91II92I 



|9i + 



^| (*i + • $>{k 2 ){q 2 x A d (g 1 ))(9 _ 2 x A e (g 2 ))|| . 



(30) 
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(2 2) 

Proceeding analogously for FqI we obtain 



F ™ = s t * * + *> U£S-a 1 M ' (31) 



A direct computation of this object turns out to be extremely cumbersome. We will need 
to wait until sec 



IV 



('2 2) ( 2 2) 

where we will be able to relate Fq£ with a known term of Fqj . Its 



explicit expression in terms of the A fields can be found in Eq. (84). 

We have thus obtained the wave functional to 0(e 2 ) by extending the method first devised 
in Ref. pQ to the next order. The different contributions to ^gl^] are summarized in Eqs. 



(15), (22), (30) and (84). This result satisfies the Schroedinger equation by construction. 
It is also explicitly real. On the other hand, we can not claim (a priori) that the Gauss 
law is satisfied, as it has only been used in some intermediate computations. At 0(e) it is 
possible to directly check that the Gauss law is satisfied. A direct check at 0(e 2 ) turns out 
to be extremely difficult to obtain, due to the complexity of the expressions involved. In sec. 



IV we will devise a method to test the gauge invariance of the expression obtained in this 
section. Finally we want to stress that the computation we have performed in this section 
has been carried out without any regularization. The final result happens to be finite but 
formal manipulations have been performed on potentially divergent expressions. We will 
come back to this issue in sec. |IVj 



III. DETERMINATION OF ^ GI [A] 

In the previous section we have been able to compute the ground-state wave functional 

at weak coupling at C(e 2 ). Nevertheless, the procedure is difficult to automatize. First, 

regularization issues have been completely skipped in the previous computation and, second, 

the Gauss law is implemented in a partial, and somewhat ad hoc, manner. This last problem 

could be overcome by reformulating the Schroedinger equation in terms of gauge invariant 

variables. One such formulation was originally worked out in Refs. j5H9j^] Here we mainly 

follow Ref. [2j, where a modified approximation scheme was devised. Even though the 

original motivation of those works was to understand the strong coupling limit, it is not 

5 In those references the regularization of the Schroedinger equation was also addressed, dealing then with 
the other potential problem of the computation of sec. |XXJ. 
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difficult to see that the approximation scheme worked out in Ref. [2] could be reformulated 
to provide with a systematic expansion of the weak coupling limit. We do so here and 
compute the ground-state wave functional to 0(e 2 ). The only relevant information for us 
will be the change of field variables used. The initial new field variables will be complex: 

A := l(A 1 + iA 2 ), A:=±(A 1 -iA 2 ) . (32) 

Therefore it is convenient to also change to complex variables the space and momentum 
components in the following way (note that k and z are defined with different signs): 

z = x\ — ix 2 , z = X\ + 1X2, 

k = -(k\+ik 2 ), k=-(ki—ik 2 ), k-x = kz + kz, (33) 

d = ±(d 1 + id 2 ),8=±(d 1 -id 2 ), dd=- A W 2 . 

A and A are still gauge-dependent degrees of freedom. Instead, in Refs. [SHH] new gauge 
invariant fields named J were introduced. We will then use the following change of variables: 
(Ai,A 2 ) — > {A, A) — > ( J (A, A), A(A, A)) where the relation between both variables is the 
following: 

A a = A a (34) 
J a = 2i (M^) ac A c + 2 - ({dM^)M^ l ) a , 

where M< is an invertible matrix, which is a function of A, defined implicitly by the equation 

A = ^M t - 1 <9M t , (35) 

which inverted yields (for a more compact expression see Eq. (5) of [7]) 

Mix) = 1 - e4-(BA) + e 2 ^-3A^-dA + C(e 3 ) , (36) 
V 2 V 2 V 2 



M\x) = 1 + e4-(dA) + e 2 ^d ( -^dA] A + 0(e 3 

v 2 v 2 Vv 2 / 



(37) 
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These equalities naturally lead to consider the following Green functions: 



G{z) = w 5W(z) 

Oz 

G(z) = 



d 2 k 1 z 



(27r) 2 k it zz + e 2 ' 
d 2 k _ e &g\ _ 1 ; 



(27r) 2 k vr zz + e 2 



Also, a useful relation reads 



which can easily be checked in momentum space. We also need (T F = 1/2) 

(M f ) ac = —Tr[T a M^T c M^ 1 } . 
Tf 

The Gauss law operator can be written in a compact form in terms of A and J: 



I a (x) = (D ■ E) a {x) = -i I [ D^^M + Df 5JC(y) N " 



<L4 6 (x) * <L4 fe (a;) / 5J c (y) x 5A b {x) 



Not surprisingly the dependence on J drops out, since it is possible to prove that 



+ D 



8A\{x) * 5A\{x) 



= 0, 



where we have used the following properties (keep in mind that M a } = M ca ) 



5A b (x) 
5A b (x) 
SA b (x) 

de 



2iMl(y)6(y - x) , 
2 



i 8 -^MA d {y) + ((dM\y))M^(y)) r _ 



5A b (x) e5A b (x) 
e(i) (-fe bh )Ml:\x) = e(V\ f edh Mt:\y), 

\ / yx \ / yx 

G{y-x) [M^\y)M\x)] 



yx 



Therefore we obtain 



I de 



SA b (x) 
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for the Gauss law operator. 

In Ref. [9] it was already shown that it was possible to write the Hamiltonian as a pure 
function of J up to terms proportional to the Gauss law, which vanish when applied to 
physical (gauge-invariant) states. If we drop those terms the Hamiltonian readsp' 



'ir(z — w) 5J a {w) 5Jb(z) 2n J SJ a {z) ' 

which we split into H = + Hi, where is the first line and Hj the second. It 



is important to note that the last term in Eq. (49) only appears after regularization of a 
divergent integral. 

We can now obtain the vacuum wave functional in powers of e. We write 

V GI [J] = exp(-F G i[J\), (50) 

where (following the notation of [2j) 

-2F GI [J] = Jfi% 2 (x l ,x 2 ) J*{x x )r*{x 2 ) + lti% a3 (xi,x 2 ,x 3 ) J ai (^i)^ 2 (^)J a3 (x 3 ) 
+ j fS^^i,x 2 ,x s ,x A ) J a ^x 1 )J a ^x 2 )J^(x 3 )J a ^x i ) + ... (51) 

and the kernels fafa 2 (xi,x 2 ), fafa 2 a 3 (xi,x 2 ,x 3 ), etc., have the expansions 

/Safri.a*) = /o 2 i 1(l2 (xi,x 2 ) + e 2 / 2 (2 i ia2 (xi,x 2 ) + • • • 

fafa 2 a 3 ( X l, X 2, X 3 ) = fo 3 l ia2 a 3 (xi,X2, X 3 ) + e 2 /fi^^l, X 2 , X 3 ) + . . . (52) 
fa%a 3 a 4 ( X l, X 2,X 3 ,X^ = f^l^aa^U X 2 , X 3 , X4) + . . . . 



Acting with the Hamiltonian of Eq. (49) onto this expansion of the wave functional and 



equating terms of equal numbers of J's we obtain recursion relations for the kernels. These 



6 Note that in Ref. [5] the normalization of J is different. 
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read 



2— /il(xi,x 2 )+4 f /S(x 1 ,x)(fl U(x,y)/ J £(y,x 2 ) + V a6 



+e 



2vr 

2 



(53) 



6 / flfa 2 ab( X ^ x 2,X,y)(^°)ab(x,y) +3 / /^ 6 (xi , X, y) (fl 1 )^ (x, y, X 2 

for the term with 2 J's, while for the terms with p > 3 J's the recursion relation is 



-pfSU + £n( P + 2 - «)/i?.^-x.(n°)a6/£X a) 



27T 
p-1 



n=2 



+ £ n(p + 1 - n)/i B l aw _ ia (n 1 ) a ta p /£^ 1 



+e 



n=2 
2 



(p + l)(p + 2) (p+2) = Kp+1) f(p +i) 

ab\ iL lab "I „ 



aba-n 



(54) 



In these equations, we have used the abbreviations (following [2]) 



(tt°) ab (x,y) = 5 ab d y G(x,y) , 
(^) a6c (x, V,z) = ~ f ahc [8{z -y) + 8{z - x)\ G(x, y) 



V ab (x,y) = 5 ab J d z 8(z - x) d z 5(z - y) . (55) 
These equations are the same as the ones in Ref. [2] (which we have checked explicitly). Note 



that the splitting into and Hj was different there, since the last term in Eq. (49) was 
included in H^°\ 



If one were able to solve the set of Eqs. (53 54 ) exactly, one would obtain the exact vacuum 



functional, without any truncation. Therefore, those equations are a perfect playground on 
which to try different resummation schemes (as it was done in Ref. [2]). Here we focus on 
the weak coupling expansion and solve those equations iteratively. 



At the lowest (zeroth) order in e, we have to solve Eq. (53) for /| 2 ' x 2 ) with 



e = 0. Note that this equation is quadratic in f^ 2 \ thus it has two solutions. We take the 
normalizable one, compatible with perturbation theory: 



£=(509(5?! -a? 2 ) 
-V 2 



(56) 
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At higher orders it is better to work in momentum space. We define 



fi% 2 a [i (x u x 2 ,x 3 ) = exp [iS^hxi ) fi% a3 (k 1 ,k 2} k 3 ), (57 
JVvhh V i J 

/ exp [i V k iXi f ( a tl 2 a 3 a 4 ( k i, fa, fa, fa). (58 
Jh/v-h/4 \ ■ / 



The recursive solution of equations (53 54) to order e gives the following lowest order ex- 



pressions for the cubic and quartic kernels: 



f^l^aaifa, ^2, ^3) 



1(1213 



24 



(2 7 r) 2 5(fc 1 + fc 2 + fe 3 ) g^{fa,fa,fa) 



(59) 



/d 4 l 1 a 2 ;6 1 fe 2 (^i^2;gi,g2) 

where 

g (3) (fa,fa,fa) 

g^\fa,k 2 ;qi,q 2 ) 



aia 2 c -fbib 2 c 



ja\a 2 c j 



64 



(2vr) 2 5(A; 1 + A; 2 + g a + g 2 ) g^ 4 \fa, fa; q u q 2 ) , (60) 



16 (fafa(fa — fa) 



Ek ± + Ek 2 + Ek 3 \ -Efci Ek 2 
1 



+ cycl. perm. 



(61) 



Ek ± + E k2 + + E q2 

g^ 3 \fa, fa, —fa - k 2 ) =^-=^ g [ • <n- -<ii - <u) 

fa + fa 



(2fa + fa) fa (2k 2 + fa) fa 



(62) 



Eh 



- g (3) (fa,fa,-fa - fa 



Ek 2 
4 



5 (3) (?i^2,-gi -92) 



fcl + &2 

(2gi + q 2 ) gi (2g 2 + <Zi) g 2 



1-1-1 



qi + q2 
and E q — \q\. 

Note that the various are not fixed completely, since they are multiplied by local 
sources. Therefore, only the completely symmetric combination is determined, any anti- 
symmetric term would give zero when multiplied by the sources, as they form a completely 
symmetric function. 

Using the expressions for f{f \ in Eq. (53), the order e 2 -term in is given by 



A 2 l ia2 (k) = 6 aia2 ^— 2 [1 + N] 



C A k 2 



(63) 
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where 

N= w{I^i 3,3>ik ' p ' - p ~ k) + Sw.l gWik - p - ~ k - ~ p) ) ■ (64) 



It is possible to perform this integration, albeit numerically. The potentially divergent 
terms vanish after doing the integration over the phase of the complex number. We obtain 

N = 0.025999 (8rr) . (65) 

Note that it is real. This is not trivial to predict a priori since g( 3 )/( 4 ) are complex functions. 
As we will see this is a strong check of the computation. The kernels f( n \ n > 5, become 
nontrivial only at higher orders. 



Note that the results above are nothing but Taylor expansions of the analogous set of Eqs. 
in Ref. j2] to the appropriate order. In practice this means setting m = in their computation 



and adding the first term in Eq. (63). This last term will play a very important role in the 
comparison with the results of the previous section. 



Once we have an (approximated) expression for \I/g/[J] we can transform it back to the 
original A variables: \I/gj[J(.A)] = ^g/[^4]- in principle, as it is a gauge invariant quantity, it 
should be possible to write it in terms of the gauge covariant quantities B and D. However, 
since we work order by order in e, we do not need this. On the other hand, rotational 0(2) 
symmetry is preserved explicitly. 

We will use the following relation to transform J fields into A fields (where the derivatives 
are in the adjoint representation: DB = dB + e[A, B]; and we have defined J = J a T a ): 

B n J = -iM\D n ~ 1 B)M ] - 1 , (66) 



as well as Eqs. (36) and (37). 
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A. Order e 



o 



In this way at O(e ) we obtain 



-2F®[A] = - 



(k x A a (k))(k x A a (-k)) 



which is the expected free-field expression. 



B. Order e 



(67) 



At 0(e) we obtain 



F$[A\=if 



abc 



\h\k\ 



(k 1 xA a (k l ))(A b (k 2 )xA c (k 3 )) 
+ i\ (h x A a (k x ))(k 2 x i b (fc 2 ))(4 x A c (k 3 )) 
(k\ ■ A a (h))(k 2 x A\k 2 ))(h x A c (4))l • (68) 



1 / ki x k 2 + %k\ ■ k 2 



\Hk\ \(N + \k 2 \ + \k 3 \)\k 2 



This term comes from a combination of and terms, as we have to remember that 
J has an expansion in e itself. Using the invariance of the integrals under interchange of 
integration variables and the fact that the delta function allows to write one momentum 
in terms of the other two, it is possible, however tedious and nontrivial, to show that the 



imaginary term of Eq. ( 68 ) vanishes and that the real part is equal to Eq. ( 22 ) . 



C. Order e 2 

At C(e 2 ) we obtain 



-2F^ = ^f -L(* x A a (k))(k x A\-k))[l + N] . (69) 
This term is associated with the f 2 2 ^ term. 
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For the term with four gauge fields we obtain 



-2Re4 2 ; 4) 



fai a 2 c fbl b 2 c f ffa + g + - + ^ 1 x l a2 (fc 2 )) x 



1 

4 " Jk/iM,<fi,ty 1 ' 1*1 + £ 

+f a ia2 c fblb2 c f j ( £ + £ + gJ + ( g ) _L 

k 2 \\ki + k 2 \ \ki\ 
(h x l ai (*i)(*2 • l a2 (4))(^ 1 (?i) x i 62 ^)) 



r- c f lb2C / ^ { ( - yir ) i(*i x ^ ai (*i))(*3 x ^(* 3 )) 

JftMMM \i=i / ^ \l*i + *2| |* 3 |/ A; 2 £; 4 

(fc 2 • l a2 (fc 2 ))(£ 4 • - (*2 x l a2 (fc 2 ))(*4 x I b2 (fc 4 )) 

+ - - 1 - - (fci x l ai (*i))(* 2 x l a2 (£ 2 )) 
|* 2 |(*3 + * 4 ) 2 *f 

(fc 3 • ^(£3)) (4 + fc 4 ) • ^(h) - (k 3 X ^(£ 3 ))(*3 + h) X ^ 2 (£ 4 ) 

+ r ia2C f bib2C [ t(i2*i\ — - h^Jk^ _ 

2 1 



■|*3 + *4||*i| (£; 3 + *4) 5 

(h X l ai (*l))(*2 X l a2 (fc 2 ))(i bl (^3) X i b2 (fc 4 )) 



_ 2fai a 2 c f b lb2 c I £ If^kA (h X l ai (*l))(*2 X A^(k 2 ))(k 3 X ^(£3)) 

;£ 2 x ^(h) 



'Vl:V2.V3,V4 \i=l 

1 / 1 



(N + |* 2 | + |A; 3 + h\)\h\ \\h + * 4 |*1 

I— =(k 2 x (* 3 - h)(h ■ ^(h)) + h ■ (k 3 - h)(h x ^(h))) 



+ 



|A; 3 + fc 4 |fc2fc, 

- r \- 2 (h x * 3 (* 4 • ^ 2 (* 4 )) - h ■ h(h x ^ 2 (* 4 ))) 

|*2|*3* 4 

+ ^ , ~ 1 - , ~„ (*1 X fc 3 (*3 + * 4 ) • ^ 2 (* 4 ) - *1 • * 3 (*3 + * 4 ) X ^(£4)) 
1*211*3 + * 4 1 2 H V 

4 



fai a 2 c fblb2 c f $ \J2ki) (*1 X A ai (*l))(*2 X l a2 (^ 2 ))(4 X ^(£ 3 ))(£ 4 X ^(h)) 



'4 V i=l 



(Ei N)(N + 1* 2 | + 1*3 + * 4 |)(N + 1* 4 | + l*i + *2|)l*ilN 

\k\q{ - (fci X fc 2 )(gi X cf 2 ) 



|* 2 ||* 4 |(*l + *2) 2 
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N L ( 2 Qi-Q2 + x fea)(gi x g 2 ) \ L _ N + N + \h + k 2 \ \ 

h + k 2 \\ V <?1 J k\? 2 J\ |fci + fc 2 | / 

/ ( 2 ki-k 2 \ ( qx -q 2 \_ ^ {k x x fc 2 )(gj x g 2 ) \ / _ ^ + fej + \ki + 

vv § A ^ ' A 1*1+^1 



-2iImF ( § 4) = 

! 

paia2C fb\\> 2 c 



if'^'-f""^ [ $ ( V/.- ) (/. : , x .ra : :!!i4 x ^(^))(^(4) x 

\ i=1 J 



k( + 2fci • A; 2 fci + fci • k 2 \ 1 



(Ifc'il + l^l + lAii + AiaDlAiillfcal V \ k i + k 2 \\k\\ \k\ + k 2 \ 2 J k 2 \ \k\ + k 2 \ \k\\ 
+tf a ia2 c fblb2 c f j \J2ki) (h x l ai (£i))(^2 x A a *(k 2 ))(k 3 x ■ ^(h)) 

JvlMmm \ i=1 J 

1 ( 2k\k\\k x + k 2 \ - \h\kKh + k 2 f - kflhKh + k 2 f + kf\k 2 \(ki + k 2 f 



k-^ k<2 /c^ (y A^3 ~\~ A^)^ 

+k 2 1 \k 2 \h-(k 1 + k 2 )) 
+2 1 f fcg • (2fci + k 2 ) h-h | fci-fca 



(N + N + I^s + ^DN I |A; 3 + A; 4 |A; 2 A; 2 |A; 2 |A; 3 2 £; 4 2 |fc 2 | |fc 3 + A; 4 | 2 A; 2 

+iri a 2 c fblb2 c I $ (i^ki) (h X A^(h))(k 2 X A a *(k 2 ))(k 3 X ^ (fc 3 ))(fc 3 • ^(h)) 



1 1 1 - - 

- - ^ - - 2^ ^ — — — — — ^——^ — —(ki ■ k 3 ) 

\k 2 \(k 3 + k 4 )*kl (\h\ + \k 2 \ + \k 3 + fc 4 |)|fci| Nl** + ktfk 2 

+ifai a 2 c fblb2 c f $ ( Vfci ) (fci X l ai (^l))(^2 X A a2 (/T 2 ))(£ 3 • ^(hm, X ^ 2 (£ 4 )) 



1 



4 



+2if a ia2 c fblb2 c I nf^kAikiX A^(ki))(k 2 X l a2 (fc 2 ))(4 X X ^ 2 (£ 4 )) 



1 - 

X fc 3 ) 



(|fci| + N + |fc 3 + fc 4 |)M Ml** + A; 4 | 2 fc 3 2 

+2ifai a 2 c fblb2 c I $ [Y^kA (ki X A^(ki))(k 2 X l a2 (4))(4 X ^(4))(^4 X ^ 2 (£ 4 )) 
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2(£x x k 2 )^\k 2 \\k 3 \{k, 3 ■ h) - \h\lh\ 3 + (h + k 2 f (2k 3 ■ k 4 + ktj | 

(Ei \h\){\ki\ + N + \h + h\)(\k 3 \ + \h\ + \h + hl^kilkilhlkjih + k 2 ) 

2(h x k 2 ) (2k 3 ■ k A + £f) 

{\h\ + \k 2 \ + \k 3 + h\)(\k 3 \ + \k 4 \ + \k t + ^Dl^il^l^l^ll^i + h\ 
If — 2fc 2 x fci fci x fc 3 fci x A; 3 



(|^i| + \h\ + \k 3 + hDlh] { \k 3 + h\kjkl \k 2 \kjkl \k 2 \\k 3 + k 4 \ 2 kj 



(71) 



The last two equations can be rewritten in several ways, yet, without an organizing 
principle, their sizes remain more or less the same. 

The resulting expression for the ground state wave functional seems to have a non- 
vanishing imaginary term. This is at odds with expectations, and with the result of the 
previous section. The real part does not look at all as the result obtained in the previous 
section either. We discuss this puzzling situation in the next section. 

IV. COMPARISON BETWEEN BOTH APPROACHES 

If we compare the expressions we have found for the ground state wave functional in 
Sees. [IT] and III we see that they look completely different. Even more so, whereas ^gl is 



explicitly real, ^gi has, a priori, a non-vanishing imaginary term. Only the O(e ) expressions 
are trivially equal. Starting at 0(e) we can get agreement between both expressions after 
quite lengthy and non-trivial rearrangements. 

At £>(e 2 ) a direct comparison by brute force turns out to be completely impossible. In 
order to compare expressions we need an organizing principle that allows to split the com- 
parison into pieces. The procedure we follow is to rewrite ^gl hi terms of J and A (actually 
we will use the variable 6 defined belowQ). If ^>gl and ^>gi are going to be equal, all terms 
proportional to A (or 9) should vanish. Moreover, to a given order in e the polynomial in A 
is finite so only a finite number of terms need to be compared. 

In order to perform such comparison to 0(e 2 ) we need the following relations: 

M f = e ed = 1 + e6 + j6 2 + C(e 3 ) , (72) 



7 The field 9 could be interpreted as a kind of generator of complex SL(N, C) gauge transformations, see 
Ref. 0. 



21 



.4 



A a (k) 



-e9+-9 2 + 0(e 3 ), 

-M^JM* + -M^~ 1 8M^ 
2 e 

J-e[9,J] + j[9, [9, J]] 



+ d O-l[O,d0\ + ^[9, [9, 



£>(e 3 



J a (k) + ik9 a (k) + -f abc / 6\k - q)J c (q) - -f abc j q6\k - q)9 c (q) 



A = -M ] ~ l dM ] = 89- ^[9,89] + ^ [9, [9, 

6 — o. 



0{e 3 



1 — ' r ^- le cabc / nb/t ^\ TCI -h le eabc 



7 P 

| jbcd jdea 

+0(e 3 ) , 



bed rdea 



0\k - «f- p)J c {q)9*(p) - — f bcd f 



9 b (k-q 



A a (k) = ik9 a {k) - ^/ a6c / gfl fe (£ - q)9 c (q) - ^ f bcd f dea / [Jfeg - kq] 9 b (k - <f - v)9 c (q)9 e (p) 



(73) 

(74) 
(75) 

: (g)r(p) 
(76) 



+0{e 3 



(77) 



where = -i9 a T a . 



For the (9(e°) and the 0(e) contributions of Fql it is possible to show that the 9 terms 



vanish and the rest agrees with Fqj in a direct fashion by just inserting the relations (76) 
and (77) into FqI and FqI and summing coefficients of terms with equal numbers of J's 



and #'s. However, for the 0(e 2 ) contributions, even after these simplifications, a brute force 
attack on the problem leads to expressions too large and complicated to directly show the 
equivalence between both expressions. 

At this respect it is better to use some intermediate expressions of the ^gl computation 
that better agree with the structure of the ^qi result in terms of J. Particularly relevant for 



us is Eq. (|29j), which relates F G \ 4) with (6F%l)/(6A). We can write F G 1] L [J, 9] = F^ L [A(J, 



in terms of g( 3 \ Using 

5A b {q) 5 fSA b {q) 



5A?(p) J q 5Am 5A b (q) J q SAf{f) 6A b (q) 

f SA\q{) 5J c (q 2 ) 5 f SAbfa) ( 5J c (q 2 ) 5 5 \ 

SA i(P) 5Ah &) 6JC &) h* 6A i(P) \8A b {&)H c (&) q2) 5A b (q 2 )J 
= I Vu + iS 2i ) («)^ + I (S U - iS 2i ) (-«|^ + ^) + 0(e) , (78) 
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we have 



8At(p) 



fx , -x s^fxM^ ( p g( 3) {k u k 2l p) k\ , . 
[Su + i5 2 i) + ^ ~ I ~p ^ + 2~^T I r - ; ''-'J ) J "<*'-> 

+ + « ^ - - - jfo) + JL + fc 2 ) + k 2 ) J ^(^(Afc) 

+|(tf u + i$H)2^ - (<J M - i^)2^|^(^i)^ 2 (^)} + 0(e) . (79) 

With this we can write Fq£ [J, 9] as a second order polynomial in g( 3 \ This gives us the 
guiding principle to try to reconstruct which is also a second order polynomial in g( 3 \ 
This term should be proportional to J 4 and we find that indeed it is. 



In Eq. (29) one can see that all terms in Fx,} '\J, 9] have a prefactor of , r , , r * — , — T . As we 



need the gauge (9) dependent terms to cancel with the corresponding terms from F G °l and 
FqI, that don't have this prefactor, we find a second guiding principle, which is to rewrite 
the 9 dependent terms of F G 2 ^ [J, 9} in such a way, that this prefactor drops out and then try 
to find a form similar to the gauge dependent contributions of FqI and Fq\. To do so we 
extensively use the Jacobi identity and the invariance of the integrals under interchange of 
integration variables, as well as the delta function. We also use the fact that the integration 
kernels can be taken to be completely symmetric under the interchange of the variables of 
two equal fields (for instance J ai (k\) J 0,2 (k 2 )) . Still the computation is very lengthy and we 
will give some details in a different publication. In the end we obtain 

4°2 = \ [ tj\k)J\-k) + e [ tlj^ki) ^f^J a (h)9 h (k 2 )J c (h) 
1 H \k\ Jk/i,k/ 2 M \~l J \k 3 \ 

1 J*fi,*fi,<$.,qi \ i J \ \ki + k 2 \ \ki\ J 

+e 2 f a ia2Cfblb2 c I f (^k, + ) J^{k l )9^{k^{q l )9\q 2 ) 
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x ( \^r-fa(fa<l2 ~ faq 2 ) + p ] ^M i + 92X9291 - 9 2 9i) 

I9l + 92| 



. 2e 2 /ai a 2C/6l 6 2C / £ ( + gl) j 9ai i)d a 2 2)9bl ^ i)9b2 ^ 2) ^m2qi + O(e 3 )(g0) 

J*/l,ti2,<fL,(f2 \ i J \fa+fa\ 

-r bc I *(£, A g(3)( ^ ,fa ^i)Afo^3) 

-y abc I $ if. A k ^nk 1 )6\k 2 )6\h) 

JftMM \ i=1 J \fa\ 

f (£& + «)) g(3)(fcl '^~ fcl - ^W i)^) Aft)^) 



a\a 2 C 



-ef aia2C f 



2 



Jth.,*k,<jL,(£ \ i= l J 

q 2 qi (fa, fa, -fa - k 2 ) q 2 k\ 



x 



2 

(a\a 2 c f b \\> 2 c 



_^ja\a 2 c jb 



91 + 92 16 (gi + g 2 ) 2|fc 2 

/ t + f - I-) ^(^1)^(^)^(91)^(92) 

J^Mi,^,qti \ i= i J \\fa + k 2 \ \kx\J 



Jyu%M,<t \ i= i J 

( fa It - J x , , (fa + fa) _ 

x ^(^92 - faq 2 ) + 4^ p-9i92 

\\fa\ \fa + fa\ 

ai a 2 c fblb2 c f i I ST^(t , -*\ \ na U t \oa 2( t \ nh f ^\ nh 2 f ^\ fafa<M2 



+4ef aia2C f 



$ \J2(fa + ft) I ^ ai (^i)^ 2 (^)^ 1 (9"'i)^ 2 (9"2)^p^| + 0(e 2 )(81) 



-^f aia2C f lb2C I $ (j2(fa + ti)) 9 {4) (fa,fa; 9i,92) J ai (fa)J a2 (fa)J bl (Qi)J b2 (Q2) 

+ y ai a 2 c fblb2C f f f E(£ + 9i) ) 9 (3) (fa,fa,-fa - fa)J ai (fa)J a2 (k2)J bl (qi)6 b2 (q2) 

I $ [jZ(fa + ®) f i ^ ± rr - ^) ^(^1)^(^2)^(9-1)^(9-2) 

J^iM,^,gfe V i / \\fa + fa\ \fa\J 
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aiQ2C rbib 2 c 



yaia2C j 



t T,(ki + ft) J ax {k)J a2 {k2)e H {qi)o h \Q2) 

q x q 2 g ( - 3 \h,k 2 ,-ki - k 2 ) 



91 + ft 
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_j_2 y«i«2C yfei&2C 
2 y a l a 2C jb\b 2 c 



/ ^ + ft) ftft 



ft ^2 
ft + 92 2|jfe 2 | 

fci + A; 2 fci 
|/ci + fc 2 | \ki\ 



/ ai (£i)0 O2 (4)0 6 Hft^ 62 (ft) 



I I feft + fi) ) ^ 2gl i' A" 1 (fei)^ a2 fe)^ 1 (ft)^ (ft) + Q(e) ■ (82) 



Ifcl + fc 2 



(2 2) 

We now move to , which is associated to a one-loop computation. We have already 
mentioned in sec. [il] that its direct determination in terms of A fields is not feasible. Again 

(2 2) 

we follow the strategy of rewriting F G £ in terms of J and 9. Actually what we use is Eq. 



(82), which we plug into Eq. (31) after we rewrite the functional derivatives in terms of J 



and A using Eq. (78). The calculation simplifies a lot and we find 



7(2,2) 

GL 



CA 
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J_ (±g®(k, Pl -k-p) + \ 1 9 (4) {P, fc; -p, -k) ) J a (k)J a (-k) + 0(e) . 
\k\ \P 2p 



*3) 



(2 2) 

This result allows us to write F G l in terms of the gauge fields. It reads 



F% 2) = —N— [ -i-(k x A a (k))(k x A a (-k)) 



GL 



4vr \k\ 



(84) 



where N has been defined in Eq. (64) 



We can now combine all the different contributions (in an, again, not completely trivial 
computation). We obtain the following equalities 



F G l[A(J, 



n 2 r u2 

F G i[J] - ^f- / -=r- J a (k) J\-k) + 0(e 3 ) 



4tt Jfi \k\ 2 



55) 



or in terms of the gauge fields 



F GI [J(A)} = F GL [A] - [ * (£ x A"(k))(k x A\-k)) + 0(e 3 ) 

4vr J* \k\ 2 



(86) 



25 



The first equality implies that F G l[A] is gauge invariant to 0(e 2 ), the second that F GI [J] is 
real to 0(e 2 ). We stress that FqI, FqI, and Fq^ are real, which is not evident at all as 



written in Eqs. (|80j), (|81| and (|82j). 

Overall we get complete agreement except for one bilinear real extra term in Fqi- Its 
origin can be traced back to the appearance of the last term of the Schroedinger equation 



in Eq. (49). In turn this term appears from an anomaly-like computation only after the 
kinetic operator has been regularized. Note that Fql was obtained without regularizing the 
theory, working with formal expressions. Therefore, even if the final result was finite, one 
could have missed contributions of this kind. For the other terms of F we have got a double 
check, which gives us strong confidence in our result. 



V. CONCLUSIONS 

We have computed the Yang-Mills vacuum wave functional in three dimensions at weak 
coupling with 0(e 2 ) precision. We have used two different methods to solve the Schroedinger 
functional equation: (A) One of them generalizes to 0(e 2 ) the method followed by Hatfield 
at 0(e) pp. We have named the result obtained ^^[A]. (B) The other uses the weak 
coupling version of the gauge invariant formulation of the Schroedinger equation and the 
ground-state wave functional followed by Karabali, Nair, and Yelnikov [2]. We have named 
the result obtained \&g?j[.A]. Each method has its own strengths and weaknesses, and they 
are to some extent complementary. 

The computations performed with method (A) are relatively simply and the results are 
explicitly real. The generalization to four dimensions of the 0(e 2 ) computation do not 
present major conceptual problems. Note that this is the order at which we expect to 
start to see the running of the coupling constant in D = 4. On the other hand, such a 
computation has two major drawbacks. First, the implementation of the Gauss law is not 
done in a systematic way but partially implemented in some intermediate steps. Therefore, 
we cannot guarantee a priori that the final result satisfies the Gauss law. As the results 
rapidly grow in size and complexity, a direct check turns out to be unfeasible. Actually we 
were only able to check the Gauss law with the help of method (B). The major drawback 
however is that the computation has been performed with an unregulated kinetic operator. 
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Whereas all computations can formally be carried out obtaining a finite result, some terms 
may be missed in this way. 

The computations with method (B) are somewhat more involved. Rather lengthy ex- 
pressions appear when we rewrite the wave functional in terms of the gauge fields A, which, 
moreover, look complex. Trying to prove by brute force that the result is real turns out 
to be impossible. Actually we only manage to prove it after a careful comparison with the 
result of method (A). Moreover, a possible generalization to four dimensions does not look 
trivial. On the other hand, method (B) is particularly appealing, as it directly works with 
gauge-invariant degrees of freedom. Therefore the Gauss law is automatically satisfied and 



it is not necessary to impose this constraint. Note also that the set of Eqs. (53), (54) can 
be solved recursively. Therefore, it could be possible to automatize it and obtain the wave 
functionals at higher order with a combination of algebraic/numeric programing. Finally 
and most important is that the operators in this computation have been regularized. This 
produced nontrivial contributions. 

We have compared both results. It is very difficult to show that they are equal in a direct 
way. The strategy we follow helps a lot, yet it continues to be extremely complicated to show 
the equivalence of two given expressions. As we have already mentioned, this comparison has 
allowed us on the one hand to prove that ^gl was indeed gauge invariant and on the other 
hand that ^gi was real. Most interestingly the agreement between both results is almost 
complete except for one extra term that appears with the method (B). This term shows up 
from an anomaly-like computation once the theory is regularized. Such a contribution does 
not show up in method (A). Apparently, this is due to the fact that no regularization was 
used in this computation. This result is potentially very interesting because it is precisely 
this term that produces the mass gap and a linearly rising potential in the strong coupling 
limit in Ref. [9]. Therefore, it's important to understand how (and if) such a term can be 
generated in a regulated version of the Schroedinger formalism in terms of the gauge fields, 
as this contribution has not been checked with an independent method so far. However, 
as regularization in the Schroedinger formalism with gauge variables is, to a large extent, 
uncharted territory, this requires a dedicated study beyond the aim of this work. Yet, we 
plan to address this issue in the near future. 
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